Microwave spectroscopy of Josephson junctions in topological superconductors 
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We consider microwave spectroscopy of Josephon junctions composed of hybridized Majorana 
states in topological 1-D superconductors. We point out how spectroscopic features of the junction 
appear in the current phase relation under microwave irradiation. Moreover, we discuss a way 
to directly probe the nonequilibrium state associated with the 4-7T periodic Josephson effect. In 
particular, we show how the microwave driving can be used to switch from a 4n to a 2tv Shapiro 
step in the current voltage relation. 
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Josephson junctions in topological superconductors 
differ from conventional superconductor junctions in sev- 
eral fundamental ways. [IH3] Their low-energy Andreev 
bound state spectrum can generically be well-separated 
from the continuum spectrum above the superconduct- 
ing gap. Moreover, the fermion parity structure of the 
states is different, which is reflected e.g. in wave function 
overlaps. These differences are due to the bound states 
being formed through hybridization of Majorana states. 

Several distinguishing features of these systems have 
been proposed, [U [3l [5j-[9] and first experimental evidence 
of the relevant physics was recently obtained in semicon- 
ductor nanowires. [101 111] One major generic feature is 
the 4-7T periodic ac Josephson effect. [T] It is a nonequi- 
librium effect, in which the system retains memory of the 
population of the bound states during time evolution of 
the superconducting phase difference ip i— > ip + 2n. This 
results to an effective 4-7T periodic current-phase relation, 
I ((p) ~ I c sm(<p/2), pQ the consequences of which are vis- 
ible in Shapiro steps and other observables. [TJ |31 [SHI] 
Importantly, the expected double-frequency Shapiro step 
feature was recently seen in an experiment |llj . 

A well-established way to probe the spectrum of a 
quantum system (e.g. a qubit) is to drive it, and look 
for resonances as a function of the frequency of the 
drive. Driving stimulates transitions between energy 
levels and thereby also induces a nonequilibrium state 
in the system. This physics is in play in the Andreev 
bound states in Josephson junctions ^2til4] and Majo- 
rana wires. [15j[T6] Information obtained in this way can 
also be useful in characterizing the special features of 
topological superconductor junctions. For instance, a 4ir 
periodic Josephson effect is not necessarily of a topologi- 
cal origin, [17] even though it is a strong indication of it, 
but combined with additional knowledge of the spectrum 
its accidental occurrence can be excluded. 

Here, we suggest how microwave driving can be used as 
a separate control parameter for tuning the magnitude of 




FIG. 1: (Color online) Josephson junction in a ID 
(proximity-induced) topological superconductor. The two 
Majorana states rji, 772 located at each side of the junction 
are hybridized (Josephson energy Ej), whereas the remain- 
ing two 773 , 774 are spatially separated and weakly coupled 
(coupling energies Sj). The system is assumed to be current 
biased and shunted with resistance R (as shown) or phase 
biased (ip fixed), and driven with a finite frequency u)d- 

the 4-7T periodic Josephson current, which does not require 
crossing topological transitions and thereby adjusting the 
level structure of the system, which may complicate in- 
terpretation of the results. We also discuss how the level 
structure of the junction is spectroscopically reflected in 
the dc current in phase- and current-biased situations. 

Model. We describe the physics of the Joseph- 
son junction in topological superconductor (TJJ) by a 
Bogoliubov-de-Gennes Hamiltonian HBdait)- The time- 
dependence originates from the time dependence of the 
superconducting phase <p(x, t); we assume (see Fig. [I]) 
that to the left of the junction 4>(x,t) — —ip(t)/2, x < 0, 
and to the right, 4>(x, t) = (p(t) /2, x > 0. Such a Hamilto- 
nian can be conveniently rewritten in the corresponding 
instantaneous Fock eigenbasis, [TS] 
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model; H DH |2D] the results in Fi S s - H a )>( b ) show 
how Eqs. p]), Q capture the low-energy features. [28] 
The connections to the continuum at e > A are not ex- 
ponentially small as shown in Fig. [2|c) , and have a rela- 
tively weaker dependency on p. The energy gap A — Ej 
depends on the transparency of the junction. [2] 

Electromagnetic drive couples to the junction by in- 
ducing a voltage V(t) = — (sdhojd/e)sm(uidt) across the 
system, equivalent to a superconducting phase difference 
ip(t) = ^fdtV(t) = 2a d a»(u d t). 

We concentrate on the population dynamics of the 
current-carrying low-energy levels, and treat the coupling 
to the continuum as a perturbation. [7] For simplicity, 
we assume the continuum connections M have roughly 
constant order of magnitude in the relevant energy range 
around ie + fiud [cf. Fig.[2fc)], that the continuum den- 
sity of states is steplike, 7V(e) = Af9(\e\ — |A|), and that 
the quasiparticle population in the continuum is negligi- 
ble. The resulting master equation for the states 1 and 
2 is 



FIG. 2: (Color online) Spectra and connections computed 
in a Rashba nanowire model, (a) Andreev spectrum (e n > 
shown for n < 18). (b) Connections Mi,2 (black line) and 
Afi,-2 (light line) as a function of phase difference. Model 
connections from Eq. Q (dotted lines) agree closely with the 
numerical results, (c) Connections M2, n> 2 < n < 18 between 
a low-energy Andreev level and the continuum. The other 
low-energy level is weakly coupled to the continuum, |Mi jn | <C 
|Af2, n |- (d) Schematic many-body spectrum with Landau- 
Zener (Plz) and continuum transitions (Pd) indicated. Here, 
\n\n2) denote the occupation numbers of levels 1 and 2. 



where M mn (ip) 



, p) are the 



connections between the instantaneous single-particle 
eigenstates. T3 is the third Pauli matrix (charge den- 
sity operator) in the electron-hole space, and we choose 
a basis such that d_„ = d\, e_ rl = — e n . 

The low-energy physics is captured by a model Hamil- 
tonian for the hybridization of the Majorana states, [T] 

H = iEj cos (|j 771772 + iSiViVa + ifomm , (2) 

where r)j are the four Majorana operators in Fig. ^ The 
resulting connections are M^i = M 2 ,2 = 0, and 



M ±2 ,iO) 



Ej 



5% + {^) cos%/2) 



-n(f), (3) 



where 5± =61^82- Away from the level crossing at ip = 
7r, these matrix elements are proportional to Majorana 
state overlaps, which are exponentially small when the 
length of the wire segments L/2 is large compared to the 
superconducting coherence length £ = Hvp/A. 

To compare with a specific microscopic theory, we 
consider the connections in a semiconductor nanowire 



p = C{p} = -i[H' ,p\- [Me*)+A(e k y] (4) 
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H o = H o + l E [A{e k .) - A{e k )*]d\d w , (5) 
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where = HJ\f\M\ 2 siu)^ is the transition rate, P = 
( — l) d i (i i+ d 2 £i 2 j s ^he Fermion parity, and H the part 
of Eq. ([!]) involving only levels 1 and 2. KeA(e) « 
7T^r_D^(|u;d| — |e — A|). E c is a cutoff energy, originating 
from the fact that M ~ e _1 decays at energies e ^ A; 
the results below are insensitive to it. Tq is the inverse 
lifetime of the continuum levels, which are assumed to be 
better coupled to external leads than the localized low- 
energy bound states. For A + e 2 > fioJd > A — e 2 [see 
Fig. (2jd)], the result describes quasiparticles on level 2 
absorbing energy from the field and escaping to the con- 
tinuum, leading to depopulation ("cooling"). [13] The 
opposite emission process is limited by the low quasipar- 
ticle population in the continuum. 

Current-phase relation. Consider first the phase bias- 
ing condition, ip(t) — 0o + 2sd cos(uj d t), in which the DC 
part 0o is kept fixed. In this setup, in periodic nonequi- 
librium effects are not visible, but one can study spectro- 
scopic features of the junction. 

We augment Eq. Q with quasiparticle poisoning 
(parity non-conserving), relaxation (parity-conserving), 
and dephasing described by phenomenological rates, for 

10- 4 Ej/h and T d ~ 



which we assume values T q ,T 
10~ 3 Ej/h, respectively, 
rent (d v H) is shown in Fig 
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2 , S~/Ej — 10 3 . The dashed line indicates the threshold 
A — 62 ((f) for continuum transitions. The n-photon 



FIG. 3: (Color online) Time-averaged current-phase relation 
7(</>o) as a function of excitation frequency U)d, keeping AC 
amplitude fixed at hsdOJd = O.lEj and F_d = 10 _5 -Ej (lighter 
color: larger current). We take A/Ej = 1.5 and 5+/Ej = 
10 

resonances nhujd = £2 (<p) — ei (if) between the Andreev bound 
states for n — 1,2 are indicated. The frequency-dependent 
threshold near ip = n is due to LZ transitions, which play a 
role when 4>q + 2sd > n. 



contacts [13], resonant transitions can here be identified 
with sharp dips in the current-phase relation at n-photon 
resonances nfuod — 62^) — ^i(ip)- The transition rate is 
proportional to the connections M ~ 8/ Ej (forbidden [4 
if there is no overlap, 5\/ 2 — 0), but it is balanced against 
the small rates T q and r r . The second apparent feature 
is that at A + £2( < p) > h^d > A — €2(<p), the transitions to 
continuum [T3] depopulate level 2 and thereby increase 
the current. Note that in conventional quantum point 
contacts, the continuum excitation gap A — e 2 is zero at 
ip = 0, whereas here it remains finite for all ip, reflecting 
the qualitatively different energy spectrum. 

Tuning the An periodic Josephson effect. The 4-7T peri- 
odic Josephson effect [TJH] under dc bias ip(t) — 2eVt/H 
in this system requires (i) Landau-Zener (LZ) transi- 
tions at ip(t) = 7r + 2im and (ii) a large gap between 
the continuum and bound state spectra [see Fig. [2jd)], 
which allows correlations between ip and ip + 2tt be pre- 
served. The former requires a transition probability 
P LZ ~ e -8ir8 2 /(EjV) _ L 5 j The latter can be modified 

through the "cooling" effect discussed above. An addi- 
tional high-frequency signal A + Ej > liujd > A — Ej 
causes a transition to the lower-energy states with prob- 
ability Pjj > during a cycle from ip = 2nn — ir to 
ip = 2nn + 7T. 

The probability Pd can be found by considering the 
time evolution of Eq. Q in the voltage-biased case. 
Given initial condition p{t x ) = |H)(H| (or |01)(01|) at 
(p(ti) = — 7t + j3, at p(t 2 ) = 7r — j3 we have Pd = 
1 - Poi,oi(*a) ~ Pu,n(h), which is shown in Fig. ga). 
Here, j3 > 8/Ej excludes the LZ transition. The prob- 
ability behaves as Pd ss 1 — e -cHr D /eV ^ c ^ n ag ^ e 
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FIG. 4: (Color online) Tuning the An periodic Josephson 
effect, (a) Relationship between the effective parameter Pd 
and the microscopic model, for A = 1.5-Ej and V ^> Yq tT , 
S 2 /Ej. (b) Effective current-phase relation, as obtained from 
master equation Eq. Q with faujg /Ej = 1.35 (solid lines), 
and the analytical result Eq. < |10[ ) with Pd from Fig. |4|a) 
(dashed lines), (c) The Air Shapiro kink for Plz = 1 and 
Pd = 0, 0.25, 0.5 from Eq. g. For P D = the result is 
shown as a dotted line; the solid line indicates an exact result 
without the mean-field app roximation, (d) Composite result 
from Fig. Hf c) and Eq. ( 10 1 for the Shapiro steps with param- 
eters RI c J(huj r ) = 0.2, s r = RI r /(hu) r ) = 0.5, Plz = 1 and 
different Pd- The curves are offset horizontally for clarity. 



transition rate increases, in agreement with a Landau- 
Zener type argument, and is finite for hwd > A — Ej. 

In the experimentally relevant |11| Shapiro step exper- 
iment, |23) the system is fed an additional ac current 
J r sin(w r t), and the Air periodicity manifests kink 
in the dc-current-dc- voltage relation I(V) at eV = hu) r , 
whereas the first 2-7T step resides at eV — huj r /2. Here, 
we assume the bias current h(t) = I + I r sin(u; r i) + 
Id sin(wd^) contains a further high-frequency component 
LJd 3> uj r which induces transitions to the continuum. To 
find the Shapiro steps we derive a semiclassical equation 
for the phase [33J by expanding the effective Keldysh ac- 
tion of the electromagnetic circuit in Fig. [T] to second 
order in quantum fluctuations (p q , [rTHTgl |2~4"] 



S[<p] 



, r h ■ 

dt ip 
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Ij[ V cl ](t)]^{t) 
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dtdt' 



k B TS(t-t') 
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Pj[jp cl }(t,t') 
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The current Ij(t) 



I{t)) and noise Pj(t,t') = 



4 



|({j(f), j(t')}) functional can in principle be obtained 
from the master equation. Considering the part 
5 /!Ldi'6>(i - if)P J [(p 6t ]{t,if)<p'i(tl) in the last term as 
a mean field (well-defined on time scales longer than 
the relaxation/quasiparticle poisoning time), and assum- 
ing the effect of the junction on the phase dynamics is 
small {I j < lb, Pj < 1%), the saddle point becomes (c.f. 
Ref . El) 



dt<p.(t) + Ij[p.](t) + 4 2) b»](t) = hit) 



2R 



(8) 
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e 
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The correlation function of the Ohmic environment cir- 



cuit is D R (t,t') 



2ie 2 R 



6(t - t'), D R (t,t) = 0. The 



47r periodic memory effects arise solely from the corre- 
lation function Pj, as the expectation value Ij must be 
invariant under translation tp i— > tp + 2ir. 

To find analytical results, we evaluate Eq. (Jsj) in a 
simplified model retaining only the main physics (cf. 
Ref. [22J : we neglect quantum coherence, and consider 
only the average populations p + — Pll,11 ^ Pol,Q1 , p_ = 
pio.io+poo.oo Q f -j-^g ^ w0 U pp er anc j i ower many-body states 
[see Fig. [2^d)]. The time evolution is constructed from 
LZ transitions p± i-^ (1 — Plz)p± + PlzP^ at y = 
7r + 27m, and continuum relaxation p + i— > (1 — Pjj)p + , 
P- i— > p_ +P]jp + at tp = 27m. Under this approximation, 
the I j and Pj functionals can be found, [5T] and the IV 
curve is obtained from the time average of Eq. Q . 

First, continuum transitions increase the magnitude 
of the 2tt Shapiro steps, as they remove quasiparti- 
cles from the low-energy levels. This is determined 
by the effective current-phase relation Ij[tp](t). For 
slow or small deviations around the trajectory tp(t) — 
2Vt, it can be approximated by the time-local relation 
Ij((f(t)) — Trp (V, tp(t))I(tp(t)) in the periodic steady 
state p {V, tp) = p (V, tp + 2ir). Wc find (V > 0) 

(l-P LZ )sin(f)sgn(cos(f))+P LZ |sin(f)| 



It = IrP, 



D ~ 



I - (1 - P D )(1 - 2P LZ ) 



(10) 

shown in Fig. ^b). At low excitation eRI r /(hui r ) <C 
1 the resulting first Shapiro step at 2eV/(huj r ) = 1 is 
similar to the supercurrent step at V = 0, but with a 
smaller effective supercurrent I c .i = 2eRIjI r /(huj r ). [55] 
The 47T periodic features are contained in the time av- 

(2) - RPJ{2V), and the 



r( 2 ) 

erage 7} ' 



(2) 



Here, I 



it) 



part 8I^P proportional to 1% is illustrated in Fig. 4^c 
Increasing Pd > cuts off the eV = hu r resonance in 
the IV curve. 

For Plz = 1 and Pd — > the correlation function fac- 
torizes, Pj[tp](t,t') = I*shi(tp(t)/2)sm(tp(t')/2), so that 
Eq. Q can be transformed [231|5S] to the stochastic equa- 
tion^ R d t tp(t)+lsin(tp(t)/2) =I b (t)+£(t) where X and 



the thermal noise £(t) are Gaussian random variables. 
This result is also shown in Fig. [4] for comparison. Note 
that phase diffusion due to a finite temperature also sup- 
presses the Shapiro steps, which is not taken into account 
in the figures. 

(2) 

Combining the steps from Ij and Pj , we obtain 
Fig. |4]jd), which shows how the IV curve reflects the 
change in periodicity from 4tt to 2ix when Pp increases. 
This switching is tunable, and can be used to establish 
both the size of the gap A — Ej and the way the effect 
is directly related to the noncquilibrium state. Namely, 
in junctions with only 2ir periodic Josephson effect, the 
high-frequency drive has no effect on the Shapiro steps 
within our model. 

Summary. In summary, we consider microwave spec- 
troscopy of Josephson junctions formed from hybridized 
Majorana bound states. We discuss what spectroscopic 
features manifest in a phase-biased situation, and how ac 
excitation can be used to tune the nonequilibrium state 
in the 47r periodic Josephson effect. High-frequency prob- 
ing (hujd ~ 400 mK x ks) was achieved experimentally in 
Ref.[5n]for superconducting quantum point contacts, and 
we expect high-frequency manipulation of junctions with 
Majorana bound states would also be feasible. 

This work is supported by the Academy of Fin- 
land, DFG grant RE 2978/1-1 and the EU FP7 project 
SE2ND. 
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Phenomenological rates 



Quasiparticle poisoning (g), relaxation (r) and dephasing (d) can be described with the master equation Eq. (4) of 
the main text, by adding phenomenological rate terms: [22] 



c q {p} = r q J2 x(d k P){p}, (ii) 

fc=-2 

C r {p} = T r [X{d 2 d 1 ) + X{d\d 2 )]{p] , (12) 

2 



C d {p}=T d sm 2 (^) £ X(4d k ){p}, (13) 



fc=l,2 



X(A){p} = A P Ai-±[AlA,p]+. (14) 

Similar terms can also be derived by integrating out a bosonic bath. [57] The dephasing rate is proportional to 
\d v e2(ip)\ 2 , hence the phase dependence. Assuming QP poisoning time of order 1 ps and a Josephson energy of 
Ej = 100 mK, we have HT q ~ 10~ 4 Ej. With environmental impedance R ~ lfi, the relaxation rate T r is also 
expected to be of a similar magnitude. [7] In the main text, we assumed the dephasing rate is KTd = ^ ^Ej. 
However, as long as fiT^ <C huj^^Ej, dephasing affects the results here only slightly because they do not require 
quantum coherence. 



Simplified master equation 

Below, we discuss a simplified model, which describes the quasiparticle physics of the Andreev bound states on 
a semiclassical level. This enables us to obtain analytical results for the Shapiro steps, and we can compare its 
predictions to those of the more detailed microscopic model. 

First, we assume the dephasing in the system is large enough, so that quantum coherence plays no role on time 
scales of H/(eV) (i.e. between two consequent Landau-Zener transitions). Under this assumption, the density matrix 
is projected as p <— > diagp. Further, assuming that the lower and upper two many-body levels are equivalent (i.e., 
ei w 0) the system is described by the populations p + = (pn.ii + poi,oi)/2 and p_ = (pio.io + /°oo.oo)/2- This defines 
a linear projection superoperator P{p] — {p+, P~) T ■ Effect of Landau-Zener (LZ) and continuum (D) transitions 
can then be taken into account with projected propagators U(t, t') = PTeh' dt £ (* )p T of the master equation. For 
time evolution between (p(t r ) = tt — ft i— > tp(t) = it + f3 and ip(t) i— > <p(t") = 3tt — f3 with j3 ~ S/Ej -C 2ir we have: 

U LZ ee PU(t, t')P T * ^ -/ l L z Z 1 , U D ee PU(t", t)P T * -/ d D J) . (15) 



G 



As a further simplification, we assume below that the continuum relaxation occurs instantaneously when the phase 



difference crosses the points <p = 2im. The current superoperator I{p} — P-^[d v H, p]+P T in this representation is 



d v e 2 (<p) 
-cLe 2 (f) 



(16) 



This allows straightforward approximation of the correlation functions required for the effective action. 

Of importance here is the time evolution over a single period, <p H> p + 2n (assuming tp(t) is monotonically 
increasing) : 

U c {p) = U(v + 2^p)={ UdULZ ' ° <lP<n ■ (17) 

\U LZ U D , n<p<2n 

This determines the cyclic steady state vq — P{po} of the system: U c (jp)vo(ip) = Vo(p), and the transient state 
U c (ip)vi(p) = XiVi(p). Here, the decay factor for the transient state is Ax = (1 — Pd)(1 — IPlz) an d is independent 
of p. Computing Ij = Tr Ipo = l T Ii>o we find Eq. (10) of the main text. 

It is now convenient to define the projection superoperators Qo(p) — vo{(p)l T and Qi(p) = l — Qo(<p) to the steady- 
state and transient subspaces. These satisfy Qo(p)U c (p) = U c (tp)Qo(<p) — Qo(p), Qi{p)U c (p) = U c {p)Q\{p) — 
Qo{p) 2 = Qo(p), and Qi(p) 2 = Qi(p). 
We assume the drive 

<p*(t) = 0o + 2eVt/h + 2s r [cos(uj r t + <f> r ) — cos{<p r )\ . (18) 

When eV and fotjj r are not exactly commensurate, long-time averages of periodic functions of (t) can be computed 
by averaging over </> and <j) r . The frequency component lj^ inducing the continuum transitions does not directly 
couple to the low- frequency dynamics, and is not included here for simplicity. 

We can now evaluate the time-average of Eq. (9) of the main text under drive (|18[ ). First, we remark that 
Sip c '(t") PU(t, t')P T k, within our approximations (see also Ref. H7|) . Moreover, D H (t, t) — 0, so that only one 
term contributes in the expression 

^ fl (^^^0^^a(^)^(O} = ^(*^^')<(9^)(^)n^')}<5(^-^ // )• (19) 

Evaluating now the long-time average and changing integration variables, we find 

27T J 1 r2TT 



I ( j ) [^]it) = -2R -p- I d0Tr[(V)(^ + </'o)f/(e + 0o,0o)Qi(0oK(0o)Qo(0o)]G(0), (20) 
Jo 2n J 

OO r 2TT 



e ^ f M d(t>r (l-P D ) m {l-2P LZ ) 



H ^ Jo 2-7T dt(p*(t m (9 + 2nm,0,<f> r )) 



1 1 ( Sr ^ r )2 cos(^)-(l-P I5 )(l-2P LZ jcos, , 



2V 1 - (1 - P D )(1 - 2P LZ ) 4?3 sin (^)2 + [ cos( ^ ) _ (1 _p D)(1 _ 2 P iz )]2 



(21) 



The function £*($, </>o, </> r ) is defined by p Sf (t if (9, 4>q, 4> r )) = 8. 

All information about the 47r periodicity and the drive is contained in the factor G(9). For Plz — 0, the maximum 
amplitude is obtained when the drive is 2ir periodic in 9, and for Plz = 1 when it is Air periodic. Moreover, Pjj > 
suppresses correlations between cycles, and cuts them off completely at P D = 1. 



Equation (20 1 can be evaluated analytically, in the leading order in the drive amplitude s r . In the interesting case 
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FIG. 5: Time-averaged current-phase relation for a nanowire system. As Fig. 3 in the main text, but for the nanowire 
system whose spectrum and connections are displayed in Fig. 2, corresponding to parameters given in [28]. The wire was 
chosen long, so that the energy splittings 5± of the Majorana states are small, and the resonances corresponding to transitions 
between the low-lying bound states are weak. They are more pronounced in shorter wires, as illustrated in Fig. 2 of the main 
text. Resonances associated with transitions to continuum levels are not similarly suppressed. The assumed continuum inverse 
lifetime is Y /Ej = 10" 3 . 



of P LZ = 1, we obtain 



/?> = I 



(2) 
/,() 



5lf + 0(s 4 r 



r(2) _ 
£ J,0 — 



2(1 -P D ) [(4 - 7r)Pb + 2tt] Ml 



(2 - P D )H 



V 



(ei?)- 1 /- 2 ^ 



(2) 



2(1 - P D )P^s 2 (hu r ) 3 sin ( 



tt(P d - 2) 2 (eV - hto r ) 2 {eV + tua r ) 2 (i* - 2P D cos 

4(1 -P D )P D S&HUt)* COS 



v eV 



2P D + 2 cos ( 



7vhljJ r 

eV 



2eV 



tt{P d - 2)eV(eV - hu: r )(eV + huj r ) (P 2 - 2P D cos (\ v 
(1 - P D )sl(hLu r ) 2 



2P D + 2 cos i 



KhuJ r 

eV 



(22) 
(23) 

(24) 



2) 



(Pd - 2) 2 eV{eV - huj r )(eV + hu) r ) ' 

This result is plotted in Fig. 4 of the main text. Note that it applies in the limit where the junction has a small effect 
on the dynamics of the phase tp 
approximation. 



that is, Eq. (18) is valid. Divergences (e.g. V 0) indicate a breakdown of this 



Discrete continuum spectrum 



The level spacing in the continuum part of the energy spectrum in nanowires is not necessarily small, for instance 
in a proximity nanowire setup where the effective energy gap A in the nanowire is smaller than the gap Ag of the 
proximity superconductor inducing it. Our main results apply also in this case — Figs. [5] and [6] show that only limited 
qualitative changes are expected from the discreteness of the continuum spectrum. 
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FIG. 6: Probability of continuum relaxation for a nanowire system. As Fig. 4(a) in the main text, but for the nanowire system 
whose spectrum and connections displayed in Fig. 2, and shown as a function of Sd with fixed eV = O.OlEj. The offset from 
zero is due to a finite spontaneous quasiparticle poisoning and relaxation with rates r r , 9 / 'Ej = 10 -4 . The result is not sensitive 
to continuum lifetime FfT 1 . 



